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Abstract 

Describing fusion reactions between 16 0 and 154 Dy and, between lb O and 144 Sm by 
the sd— and sdf— interacting boson model, we show that heavy-ion fusion reactions 
are strongly affected by anharmonic properties of nuclear surface vibrations and nuclear 
shape, and thus provide a powerful method to study details of nuclear structure and 
dynamics. 

1. Introduction 

Extensive experimental as well as theoretical studies of heavy-ion fusion reactions at 
energies near and below the Coulomb barrier in the past decade have offered many basic 
ideas on the effects of nuclear collective excitations such as surface vibrations and rotation, 
or particle transfer on the fusion cross section [[3]]. Recent very accurate data at ANU and 
at Legnaro are providing us with the possibility to step up to a more advanced stage of 
the research. For example, one can extract detailed informations on nuclear structure 
and dynamics from the data of heavy-ion fusion reactions. For instance, the sign of the 
hexadecupole deformation can be determined by studying the fusion barrier distributon, 
as was exemplified for the 16 0 + 186 W, 154 Sm fusion reactions p. j||]. Also, Stefanini et 
al. [Q demonstrated the important role played by the double phonon excitations in the 
58 Ni + 60 Ni fusion reactions. In this contribution we address to the effects of anharmonic 
properties of nuclear surface vibrations. Our study was motivated by the discussions 
with Jack, Nand and David. When we visited Canberra last spring, they introduced us 
their puzzle about the role of octupole vibrations in the fusion reactions between 16 0 
and 208 Pb and between 16 0 and 144 Sm. The problem is the following: in the former, a 
better agreement between data and theory concerning the fusion barrier distribution is 
obtained if one includes the double phonon coupling p]. On the other hand, in the latter, a 
good agreement obtained by a single phonon coupling is destroyed if one adds the double 
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phonon coupling p], 144 Sm has low lying quartet states which seem to correspond to the 
double octupole phonon excitations^. Therefore, it is quite puzzling that the inclusion of 
the double phonon coupling spoils the agreement bewteen the data and theory, ffowever, 
there exists a noticeable anharmonicity in the energy spectrum [0. We therefore expected 
that the puzzle raised by Jack et al. can be attributed to the anharmonic properties 
of nuclear surface vibrations. This is what we are going to report in this contribution. 
Actually important anharmonic properties are concerned not with the energy spectrum, 
but rather with the transition matrices, especially reorientation effects. 

2. Coupled-channels formalism in the IBM — a case study for ie O + 154 Dy 
fusion reactions 

In this section we make a case study for the 16 0 + 15i Dy fusion reactions by using the 
U(5) limit of the sd— IBM. It is known that the U(5) limit of the IBM corresponds to 
the anharmonic vibrator in the geometrical model of Bohr and Mottelson [|J. We choose 
the IBM to reduce the number of free parameters. 

We assume the following Hamiltonian for the total system 

n 2 

H = ——V" H IBM + Vcoupi'i'', £), (1) 

where r is the coordinate of the relative motion between the projectile and target, /i is 
the reduced mass, and £ represents the internal degrees of freedom of the target nucleus. 
Hibm is the IBM Hamiltonian for the quadrupole vibration in the target nucleus, for 
which we assume the harmonic limit 


Hibm — f-dfid (2) 

Here fid is the number operator of the d bosons and is the excitation energy of the 
quadrupole vibration. In the numerical analyses we also introduce a deviation from the 
harmonic spectrum and study its effects. 

The coupling between the relative motion and the intrinsic motion of the target nucleus 
is described by V coup in Eq. ([!]). We use the linear coupling model and the no-Coriolis 
approximation pj. We then have 

V CO up(r, £) = U N (r) + Zp ^^ + -j==f(r)Q 20 (3) 

where the coupling form factor f(r) consists of the nuclear and Coulomb parts and read 

f(r) = -Rr^ + \z P Z T e 2 ^f (4) 

dr 5 r 6 


N is the boson number, the subscripts P and T refer to the projectile and target nuclei, 
respectively. The scaling of the coupling strength with y/N is introduced to ensure the 
equivalence of the IBM and the geometric model results in the large N limit JTIJ. (3 is the 








quadrupole deformation parameter. We represent the quadrupole moment operator Q 2 ^ 
by 


Q 2 /x = + sd\ + x(d t d) (2M) 


( 5 ) 


where tilde is defined as bi^ = (—When the x parameter is zero, the quadrupole 
moments of all states vanish, and one obtains the harmonic limit in the large N limit. 
Non-zero values of the x generate quadrupole moments and are responsible for the anhar- 
monicities. This can be clearly seen from the structure of the matrix elements appearing 
in the coupled-channels calculations, which read, for example, 


Hibm + V coup — Un{t ) -\- 


Z P Z T e 2 


/ 0 F(r) 0 \ 

+ ^ - vmXF(r) - l/JV)F(r) 

1 0 /i(l - l/N)F(r) 2 e d + S- ^X^r) ) 

( 6 ) 

if we truncate by the two phonon states. F(r) is ^=/(r). We introduced a parameter 5 
in order to represent the deviation of the energy spectrum from the harmonic limit. 

We now apply this formalism to the fusion reactions between 16 0 and 154 Dy treating 
16 O as inert 0- The boson number for 154 Dy is estimated to be 4 by considering 64 and 82 
as the proton and the neutron magic numbers, respectively. The data of the E2 transition 
and the intrinsic quadrupole moment lead to /3=0.24 and y=-2.36. The excitation energy 
of the d boson, e^, has been identified with the excitation energy of the first excited state 
of 154 Dy, 0.335 MeV. 

Fig.l shows the role of anharmonicity in this reaction by comparing the fusion exci¬ 
tation function and the fusion barrier distribution calculated in four different ways. The 
dotted line is the result of the one phonon coupling in the harmonic oscillator limit. The 
dott-dashed line is the result of the one phonon coupling where the reorientation term is 
included. The dashed line is the result when one takes double phonen states into account 
in the harmonic oscillator limit. The solid line is the corresponding calculations where the 
anharmonicities are taken into account. The effects of anharmonicity can be hardly seen 
in the fusion excitation function. However, the fusion barrier distribution clearly shows a 
significant change from the harmonic oscillator limit by the anharmonicity. 

Fig.2 shows the dependence of the fusion barrier distribution on the deviation of the 
energy spectrum from the harmonic oscillator limit (the upper panel) and on the boson 
number (the lower panel). The upper panel compares the fusion barrier distribution 
calculated for four different choices of the parameter 6 in Eq. (^). It shows that the fusion 
barrier distribution is insensitive to the anharmonicity concerning the energy spectrum. 
The lower panel compares the fusion barrier distribution calculated for various boson 
numbers N. For comparison, it also contains the fusion barrier distribution calculated in 
the harmonic oscillator limit in the case, where there exist one and two phonon states. 
One clearly observes a strong dependence of the fusion barrier distribution on the boson 
number, and also the effects of anharmonicity. 

Fig.3 shows how the sign of nuclear deformation, i.e. the sign of the parameter x in 
the quadrupole moment (see Eq. ([5])), affects the excitation function of the fusion cross 
section and the fusion barrier distribution. Though the former is almost insensitive to the 









sign of the deformation, the latter drastically changes with the sign of the deformation. 
This enables us to determine the shape of nuclei with anharmonic vibrational modes of 
excitation by analysing the fusion barrier distribution. 


3. Anharmonic effects of 144 Sm on the ie O + 144 Sm fusion reactions 


We now come to the question raised by Jack et al. (see Ref. [^] for details). Our 
main concern is the anharmonic effects of the octupole excitation of 144 Sm. We therefore 
extend the model in the previous section by including the / bosons, and assume that the 
IBM Hamiltonian is given by 

Hibm = Q'fO + tfhf- (7) 

where h/ is the number operator for / bosons and e/ is the excitation energy of the 
octupole vibration. Since the study in the previous section has shown that the anhar- 
monicity concerning the energy spectrum does not so much influence both the excitation 
function of the fusion cross section and the fusion barrier distribution, we assume here 
a harmonic spectrum. The coupling Hamiltonian now includes also the coupling to the 
octupole vibration. Furher, we modify the linear coupling model in the previous section 
by treating the nuclear coupling to full order, while we keep the Coulomb coupling in the 
linear coupling approximation fllOl, [R§. The Coulomb and the nuclear parts of the coupling 
Hamiltonian are now given by 


Vc(r,0 

Vjv(r,0 


Z P Z T e 2 ,34 &Q 20 


1 + - 


_l_ 3 R'p (hQ'-io ^ 


-V a 


5 r 2 a/ 47 tN 7 r 3 a/ 47 tN J ’ 

1 + exp f— (r — R 0 — Rt{@ 2 Q 2 o + fisQso)/ V 47 riv) 


-l 


( 8 ) 


f3 2 and /J 3 in Eq. (||) are the quadrupole and octupole deformation parameters, which 
are usually estimated from the electric transition probabilities using the expression /3\ = 
An(B(E\) |) 1 // 2 /3 ZpeR^. However, this formula does not hold for anharmonic vibrators. 
Therefore, we treat /3 2 and (3 3 as free parameters and look for their optimal values to 
reproduce the experimental data. Q 2 and Q 3 in Eq. (|[) are the quadrupole and the 
octupole operators in the IBM, which we take as 


Q 2 = s ] d + sd ] + y 2 (d t d) (2) + X 2 /(/ T /) (2) , 

Q 3 = sf + X3W^ + h.c., (9) 


The model parameters are determined as follows. The standard prescription for boson 
number (i.e. counting pairs of nucleons above or below the nearest shell closure) would 
give N = 6 . However, it is well known that the effective boson numbers are much 
smaller due to the Z = 64 subshell closure ||. The suggested effective numbers in the 
literature vary between N = 1 and 3. We adopted N = 2 in our calculations, since there 
are experimental signatures for the two-phonon states, but no evidence for three-phonon 
states in 144 Sm. The parameters of the IBM Hamiltonian Eq. (^) are simply determined 
from the excitation energies of the Erst 2 + and 3“ states in 144 Sm as = 1.66 MeV and 










e/ = 1.81 MeV. The nuclear potential parameters are taken from the exhaustive study of 
this reaction in Ref. || as Vq — 105.1 MeV, R 0 = 8.54 fm and a = 0.75 fm. Finally, the 
target radius is taken to be Rt = 5.56 fm. 

The results of the coupled-channels calculations are compared with the experimental 
data in Fig. 4. The upper and the lower panels in Fig. 4 show the excitation function of the 
fusion cross section and the fusion barrier distributions, respectively. The experimental 
data are taken from Ref. The dotted line is the result in the harmonic limit, where 
couplings to the quadrupole and octupole vibrations in 144 Sm are truncated at the single¬ 
phonon levels. The deformation parameters are estimated to be /? 2 =0.11 and /?3=0.21 from 
the electric transition probabilities. The dotted line reproduces the experimental data of 
both the fusion cross section and the fusion barrier distribution reasonably well, though 
the peak position of the fusion barrier distribution around E cm = 65 MeV is slightly 
shifted. As was shown in Ref. |J, the shape of the fusion barrier distribution becomes 
inconsistent with the experimental data when the double-phonon channels are included in 
the harmonic limit (the dashed line). The good agreement is recovered when one takes the 
effects of anharmonicity of the vibrational motion into account. These results are shown 
in Fig. 4 by the solid line. This calculation has been performed using the parameters, 
@2 — 0.13, p 3 = 0.23, X '2 = —3.30, y 2 / = —2.48, and y 3 =2.87, which are obtained from 
a x 2 St to the fusion cross sections. The y 2 fit gave a unique result, regardless of the 
starting values. The non-zero y values indicate the anharmonic effects in the transition 
operators. The slight change in the values of the deformation parameters from those in 
the harmonic limit results from the renormalization effects due to the extra terms in the 
operators given in Eq. (£5[). Note that the solid line agrees with the experimental data 
much better than the dotted line. 

One of the pronounced features of an anharmonic vibrator is that the excited states 
have non-zero quadrupole moments. Using the y parameters in the E2 operator, T(E2) = 
£bQ 2 , extracted from the analysis of fusion data, we can estimate the static quadrupole 
moments of various states in 144 Sm. Here, e# is the effective charge, which is determined 
from the experimental B( E2; 0^2^) value as e# = 0.16 eb. For the quadrupole moment 
of the first 2 + and 3“ states, we obtain —0.28 b and —0.70 b, respectively. The nega¬ 
tive sign of the quadrupole moment of the octupolc-phonon state is consistent with that 
suggested from the neutron pick-up reactions on 145 Sm ||T4] . 

In the case of rotational coupling, fusion barrier distributions strongly depend on 
the sign of the quadrupole deformation parameter through the reorientation term. Also, 
as mentioned in the introduction it has been reported that fusion barrier distributions 
are very sensitive to the sign of the hexadecapole deformation parameter J2|. Similarly, 
it is likely that the shape of fusion barrier distributions changes significantly when one 
inverts the sign of the quadrupole moment in a spherical target if there exists a strong 
anharmonicity. Fig. 5 shows the influence of the sign of the quadrupole moment of the 
excited states on the fusion cross section and the fusion barrier distribution. The solid 
line is the same as in Fig. 4 and corresponds to the optimal choice for the signs of the 
quadrupole moments of the first 2 + and 3“ states. The dotted and dashed lines are 
obtained by changing the sign of the y 2 and y 2 / parameters in Eq. ([)]), respectively, 
while the dot-dashed line is the result where the sign of both y 2 and y 2 / parameters are 



inverted. The change of sign of y 2 and X 2 f is equivalent to taking the opposite sign for 
the quadrupole moment of the excited states. Fig. 5 demonstrates that subbarrier fusion 
reactions are indeed sensitive to the sign of the quadrupole moment of excited states. The 
experimental data are reproduced only when the correct sign of the quadrupole moment 
are used in the coupled-channels calculations. Notice that the fusion excitation function 
is completely insensitive to the sign of the quadrupole moment of the first 2 + state, but 
strongly depends on that of the first 3“ state. In contrast, the fusion barrier distribution 
can probe the signs of the quadrupole moments of both the first 2 + and 3^ states. This 
study shows that the sign of quadrupole moments in spherical nuclei can be determined 
from subbarrier fusion reactions, especially through the barrier distribution. 

4. Summary 

We discussed the effects of anharmonicity in nuclear surface vibrations on heavy- 
ion fusion reactions by using the interacting boson model. Our analyses clearly show 
that the fusion barrier distribution is very sensitive to the anharmonicity of the nuclear 
surface vibrations, and suggest that the puzzle raised by Jack et al. concerning the 16 0 + 
144 Sm fusion reactions can be solved by considering the anharmonicities of the octupole 
surface vibrations of 144 Sm. We have also shown that one can determine nuclear shape by 
using this high sensitivity of the fusion barrier distribution to the anharmonicity through 
reorientation effetcs. 
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Figure Captions 


Fig. 1 : Effects of anharmonicity in the 16 0 + 154 Dy fusion reactions. 

Fig. 2: Dependence of the fusion barrier distribution on the energy spectrum and the 
boson number in the 16 0 + 154 Dy fusion reactions. 

Fig. 3: Dependence of the excitation function of the fusion cross section for the 16 0 + 
154 Dy reactions and of the fusion barrier distribution on the boson number. 

Fig. 4: Comparison of the experimental fusion cross section (the upper panel) and fusion 
barrier distribution (the lower panel) with the coupled-channels calculations for 16 0 + 
144 Sm reaction. The experimental data are taken from Ref. ||. The solid line shows the 
results of the present IBM model including the double-phonon states and anharmonic 
effects. The dotted and the dashed lines are the results of the single- and the double¬ 
phonon couplings in the harmonic limit, respectively. 

Fig. 5: Dependence of the fusion cross section and barrier distribution on the sign of the 
quadrupole moment of the excited states in 144 Sm. The meaning of each line is indicated 
in the inset. 
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